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ABSTRACT 

Context. Unveiling the structure of the disks around intermediate-mass pre-main- sequence stars (Herbig Ae/Be stars) is essential for 
our understanding of the star and planet formation process. In particular, models predict that in the innermost AU around the star, the 
dust disk forms a "puffed-up" inner rim, which should result in a strongly asymmetric brightness distribution for disks seen under 
intermediate inclination. 

Aims. Our aim is to constrain the sub-AU geometry of the inner disk around the Herbig Ae star R CrA and search for the predicted 
asymmetries. 

Methods. Using the VLTI/AMBER long-baseline interferometer, we obtained 24 near-infrared (H- and ^T-band) spectro- 
interferometric observations on R CrA. Observing with three telescopes in a linear array configuration, each data set samples three 
equally spaced points in the visibility function, providing direct information about the radial intensity profile. In addition, the obser- 
vations cover a wide position angle range (~ 97°), also probing the position angle dependence of the source brightness distribution. 
Results. In the derived visibility function, we detect the signatures of an extended (Gaussian FWHM ~ 25 mas) and a compact 
component (Gaussian FWHM ~ 5.8 mas), with the compact component contributing about two-thirds of the total flux (both in H- and 
K-band). The brightness distribution is highly asymmetric, as indicated by the strong closure phases (up to ~ 40°) and the detected 
position angle dependence of the visibilities and closure phases. To interpret these asymmetries, we employ various geometric as well 
as physical models, including a binary model, a skewed ring model, and a puffed-up inner rim model with a vertical or curved rim 
shape. For the binary and vertical rim model, no acceptable fits could be obtained. On the other hand, the skewed ring model and 
the curved puffed-up inner rim model allow us to simultaneously reproduce the measured visibilities and closure phases. From these 
models we derive the location of the dust sublimation radius (~ 0.4 AU), the disk inclination angle (~ 35°), and a north- southern 
disk orientation (PA~180-190°). Our curved puffed-up rim model can reasonably well reproduce the interferometric observables and 
the SED simultaneously and suggests a luminosity of ~ 29 L© and the presence of relatively large (> 1.2 jim) Silicate dust grains. 
Our study also reveals significant deviations between the measured interferometric observables and the employed puffed-up inner 
rim models, providing important constraints for future refinements of these theoretical models. Perpendicular to the disk, two bow 
shock-like structures appear in the associated reflection nebula NGC 6729, suggesting that the detected sub-AU size disk is the driving 
engine of a large-scale outflow. 

Conclusions. Detecting, for the first time, strong non-localized asymmetries in the inner regions of a Herbig Ae disk, our study 
supports the existence of a puffed-up inner rim in YSO disks. 

Key words, stars: pre-main- sequence - circumstellar matter - accretion, accretion disks - outflows - individual: RCrA - planetary 
systems: protoplanetary disks - techniques: interferometric 



1. Introduction 

For our understanding of the structure and physical processes in 
the disks around young stellar objects (YSOs), the inner-most 
disk regions are of special importance. Furthermore, it is be- 
lieved that planet formation is a direct result of the grain ag- 
gregation and growth which should take place in the inner few 
AUs of the dusty disks around these stars. As the spatial scales 
of the inner circumstellar environment were not accessible to 
imaging observations until recently, most conclusions drawn on 
the geometry of the inner disk were based on the modeling of 
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the spectral energy distribution (SED). Typically, the SED of 
intermediate-mass YSOs (Herbig Ae/Be stars) shows a charac- 
teristic infrared excess emission, which is often interpreted as 
the presence of a circumstellar accretion disk, although the 3-D 
geometry of the innermost (AU- scale) region of these disks is 
still poorly known. 

According to the current paradigm, most of the infrared ex- 
cess emission in Herbig Ae/Be stars originates from a passive 
dust disk (Adams et al. 1987), whose thermal structure can be 
approximated with a cold disk interior and a hot surface layer 
(Chiang & Goldreich 1997). To fulfill vertical hydrostatic equi- 
librium, the disk is expected to flare towards larger radii, allow- 
ing the outer disk regions to intercept more stellar light than ex- 
pected for a geometrically flat disk (Kenyon & Hartmann 1987). 
At a certain distance from the star, most often referred to as the 
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dust sublimation radius R su bh the dust temperature will exceed 
the evaporation temperature of dust (r su bi), causing the trun- 
cation of the dust disk and the formation of a dust-free inner 
hole. This scenario also gained support from the first survey- 
type near-infrared interferometric observations of Herbig Ae/Be 
stars. Millan-Gabet et al. (2001) could measure the character- 
istic size of many YSO disks, finding that the measured sizes 
scale roughly with the square-root of the luminosity L+ of the 
stellar source (Monnier & Millan-Gabet 2002). Since this is the 
expected scaling-relation for the location of the dust- sublimation 
radius (R su b\ 00 2 ), this finding supports the idea that the 
near-infrared emission is tracing mainly hot material located in 
a structure at the location of the dust sublimation radius. More 
recently, significant deviations from the size-luminosity were de- 
tected concerning, in particular, the T Tauri and Herbig Be star 
regime (Monnier et al. 2005). These deviations might be ex- 
plained either with contributions from scattered light (most im- 
portant for T Tauri stars; e.g. Pinte et al. 2008) or from a gaseous 
disk located inside of the dust sublimation radius (more impor- 
tant for Herbig Be stars; see e.g. Akeson et al. 2005; Eisner et al. 
2007; Kraus et al. 2008; Isella et al. 2008; Tannirkulam et al. 
2008). 

In 2001, Natta et al. and Dullemond et al. pointed out that the 
heating provided by the stellar radiation should considerably in- 
crease the disk scale height close to the dust sublimation radius, 
resulting in the formation of a "puffed-up" inner rim. This rim 
might cast a shadow on the more extended disk regions, possi- 
bly affecting the thermal disk structure out to hundreds of AUs. 
The shadowing effects of the puffed-up rim might already have 
been observed indirectly in the SED. Meeus et al. (2001) argued 
that the general shape of the SED of many Herbig Ae/Be stars 
can be divided into two groups, where "group I" sources show a 
pronounced 60 jim bump, while "group II" sources have a flat- 
ter SED, lacking the 60 jim excess emission. Later, Dullemond 
& Dominik (2004) pointed out that this empirical classification 
could be explained by considering both the flaring properties of 
the outer disk and the shadowing effects of a puffed-up inner rim. 
Accordingly, group I disks might exhibit a strongly flared shape, 
allowing the outer disk to step out of the shadow region, while 
group II disks are fully self-shadowed. 

Although these studies provided some important first in- 
sight, there are major uncertainties concerning both the magni- 
tude of the puffing-up effect as well as the detailed rim shape. 
For instance, it was proposed that density-dependent dust sub- 
limation effects (Isella & Natta 2005, referred to as IN05 in 
the following) as well as grain growth and dust sedimentation 
(Tannirkulam et al. 2007) could affect the shape of the inner rim, 
resulting in a curved vertical rim geometry. In a recent study, 
Kama et al. (2009) found that backwarming effects and the pres- 
ence of highly refractory grain species might result in a dust rim 
which is located significantly closer to the star than anticipated 
in earlier studies. Using radiative transfer, they show that the 
rim is not an infinitely sharp wall, but has an optically thin re- 
gion which might extend a significant fraction of the rim radius, 
resulting in a more diffuse rim morphology. Spatially resolved 
observations are required to directly measure the shape of the 
inner rim. Furthermore, various authors (e.g. Miroshnichenko 
et al. 1999) have pointed out that the infrared emission of many 
Herbig Ae/Be stars might contain significant scattered-light con- 
tributions from optically thin circumstellar envelopes or halos. 
As investigated by Vinkovic et al. (2003) and others, high- 
angular resolution imaging observations, such as presented in 
this study, provide the only method to separate the disk and enve- 



lope contributions and to solve the ambiguities inherent to pure 
SED model fits. 

One of the strongest predictions which is common to all 
puffed-up rim scenarios is the appearance of asymmetries in the 
source brightness distribution for disks seen under an intermedi- 
ate inclination angle. This "skew" in the brightness distribution 
is a direct consequence of the vertical extension of the rim above 
the disk midplane, providing perhaps the most promising way 
to observationally distinguish between scenarios with and with- 
out a puffed-up rim. In order to detect these signatures, infrared 
interferometry provides not only the required milli-arcsecond 
(mas) angular resolution, but also offers a very sensitive measure 
for deviations from point- symmetry, namely the closure phase 
(CP) relation (Jennison 1958). First closure phase measurements 
on YSOs were presented by Monnier et al. (2006). Using the 
IOTA interferometer and baseline lengths up to 38 m, these au- 
thors measured statistically significant non-zero closure phases 
on six out of 14 stars (excluding binary stars). For five of these 
six stars, the detected closure phase signals were rather small 
(<D < 5°), while for the Be star HD 45677 phases of up to ~ 27° 
could be detected. 

In order to investigate the geometry of the inner circumstellar 
environment around a Herbig Ae star in great detail and to obtain 
further evidence for or against the existence of a puffed-up inner 
rim, we have studied the Herbig Ae star R CrA using ESO's Very 
Large Telescope Interferometer (VLTI) and the AMBER beam 
combination instrument. In the following, we will first summa- 
rize some earlier studies on R CrA (Sect. 2), followed by a de- 
scription of our interferometric observations (Sect. 3) and the 
employed models (Sect. 4). Finally, we will discuss our model- 
ing results (Sect. 5) and conclude with a brief summary (Sect. 6). 



2. Earlier studies on RCrA 

R CrA is located in the Corona Australis molecular cloud and is 
the most luminous (L bo i ~ 99 L©, Bibo et al. 1992) star of the 
very young and obscured cluster known as the Coronet cluster 
(Taylor & Storey 1984). The stellar parameters of RCrA are still 
rather uncertain, ranging from a spectral type of F5 (Hillenbrand 
et al. 1992; Natta et al. 1993; Garcia Lopez et al. 2006), A5 
(Herbig & Bell 1988; Chen et al. 1997) to B8 (Bibo et al. 1992; 
Hamaguchi et al. 2005). For the distance, we adopt the com- 
monly assumed value of 130 pc (Marraco & Rydgren 1981). 

Compared to other stars in the Herbig Ae/Be class, R CrA is 
in a particularly early evolutionary phase (Malfait et al. 1998) 
and still embedded in an extended and massive natal enve- 
lope, whose emission dominates the SED from mid-infrared to 
millimeter-wavelengths (Fig. 12, bottom, left). Natta et al. (1993) 
estimated the mass of the envelope to be ~ 10 M© with an outer 
envelope radius of 0.007 pc (~ 1450 AU). In the classification 
scheme of Meeus et al. (2001), R CrA was classified as a group II 
object (Acke & van den Ancker 2004), suggesting the presence 
of a disk with a pronounced inner rim, which might self-shadow 
the outer disk. Optical polarization measurements showed a high 
degree of linear (~ 8%) and circular (~ 5%) polarization (Ward- 
Thompson et al. 1985; Clark et al. 2000), indicating scattering 
from aligned, non- spherical dust grains. The polarization map- 
ping also showed an extended (~ 10"), disk-like structure with 
north-south orientation (6 = 189 + 5°, Ward-Thompson et al. 
1985). The total mass of gas and dust in the disk was estimated 
from sub-millimeter observations by Mannings (1994), yielding 
0.02 M©. More recently, Groppi et al. (2007) used the SMA to 
derive an upper mass limit of 0.012 M©. 
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At optical wavelengths, the star is known to be highly vari- 
able both on long and short time scales (Bellingham & Rossano 
1980). Graham & Phillips (1987) reported strong variability 
in the Ha line of RCrA as well as surface brightness varia- 
tions in the nearby reflection nebula NGC 6729, possibly indi- 
cating shadowing effects caused by material in the inner cir- 
cumstellar environment of R Cr A or clumpy accretion (Graham 
1992). Variability was also observed at radio- and X-ray wave- 
lengths (Forbrich et al. 2006). The X-ray spectrum of RCrA is 
very unusual, including a very hot X-ray emission component. 
Intermediate-mass YSOs are expected to show no X-ray emis- 
sion, since they should possess neither magnetically driven coro- 
nae nor the radiation-driven winds which cause X-ray-emitting 
shock regions. Therefore, in order to explain the detected X-ray 
spectrum, Forbrich et al. (2006) suggested that the X-ray emis- 
sion does not originate from the optical/infrared source, but from 
a yet undiscovered Class I companion. The presence of a close 
companion was also proposed by Takami et al. (2003) based on 
spectro-astrometric observations which revealed a photo-center 
displacement 1 both in the blue- and red-shifted wing of the spec- 
trally resolved Ha-line. The authors pointed out that the detected 
signature cannot be explained with a stellar companion or an out- 
flow component alone, but possibly with a combination of both 
scenarios. However, as pointed out by Choi et al. (2008), it is 
unlikely that the companion proposed by Takami et al. (2003) is 
identical to the Class I companion proposed by Forbrich et al. 
(2006), since a deeply embedded Class I source would not con- 
tribute significantly to the Ha line flux at visual wavelengths. 

The accretion activity of RCrA was estimated by Garcia 
Lopez et al. (2006) from the luminosity of the Bry-line, yielding 
a relatively low accretion luminosity of L acc = 3.2 Lq, corre- 
sponding to a mass accretion rate of M acc = 10" 712 M@/yr (for a 
F5 star). Besides indications of active accretion, various outflow 
tracers have been reported for R CrA. For instance, a compact 
bipolar molecular outflow with an east- west orientation (Walker 
et al. 1984; Levreault 1988; Graham 1993) as well as several 
Herbig-Haro objects (in particular HH 104 A/B) have been as- 
sociated with RCrA (Hartigan & Graham 1987; Graham 1993). 
However, more recent studies (Anderson et al. 1997; Wang et al. 
2004), convincingly identified the source IRS 7 as driving source 
of these outflows, making a physical association with RCrA 
rather unlikely. 
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Fig. 1. wv-plane coverage obtained with our VLTI/AMBER observa- 
tions on RCrA (AT-band only). 
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Fig. 3. Closure phases recorded on 2008-06-02 on R CrA in the AT-band 
using the AMBER beam commutation device (BCD), showing the ex- 
pected CP sign change between BCD OUT (black) and BCD IN (red). 



3. Observations 

We observed R CrA during four nights in June 2008 using the 
VLTI near-infrared interferometric instrument AMBER (Petrov 
et al. 2007), which combines the light from three of the 1.8 m 
auxiliary telescopes (ATs). For all observations, AMBER's low 
spectral resolution mode (LR-HK) with a spectral resolution of 
A/AA = 35 and a wavelength coverage from 1.4 to 2.5 yum (H- 
and K-bznd) was used. The telescopes were placed on stations 
E0-G0-H0, forming a linear array configuration with baseline 
lengths of 16 (E0-G0), 32 (H0-G0), and 48 m (E0-H0). Due to 
this linear telescope arrangement and the spectral coverage pro- 
vided by the AMBER instrument, each of our AMBER observa- 
tions provides a good sampling of the wv-plane towards a certain 
position angle (PA), covering spatial frequencies between ~ 4 
and 25 x 10 6 A, allowing one to cover a wide range of the visi- 
bility function (0.8 < V < 0.1). By obtaining a large number of 
observations towards different hour angles, we also achieved a 



1 Earlier spectro-astrometric observations were presented by Bailey 
(1998), but did not show this spectro-astrometric signature. 



good PA coverage, probing PAs between 16° and 113° (Fig. 1). 
In order to adjust for changing atmospheric conditions and to test 
for possible systematic effects, we recorded data with different 
detector integration times (DITs), namely 50 ms (default value), 
100 ms, and 200 ms. 

The science observations on RCrA were interlayed with 
observations on interferometric calibrator stars, allowing us to 
monitor the instrumental transfer function. The transfer function 
is used in the course of data reduction to correct for wavelength- 
dependent instrumental or atmospheric effects. Each observation 
block (either on R CrA or on the calibrator star) typically con- 
sisted of five data sets, each containing 1000 (for DIT=50 ms) 
or 500 (for DIT=100 ms or DIT=200 ms) individual spectrally 
dispersed interferograms. 

For data reduction, we used the amdlib data reduction soft- 
ware (release 2.2). This software employs the P2VM algorithm 
(Tatulli et al. 2007) to derive wavelength-dependent visibili- 
ties and closure phases (CP). The wavelength calibration was 
done using the procedure described in Appendix A of Kraus 
et al. (2009). Following the standard AMBER data reduction 
procedure, we select the 10% of interferograms with the best 
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Fig. 2. Transfer function of the night 2008-06-09 for one spectral channel around 2.2 jim and observations with DITs of 50 ms (left) and 100 ms 
(right). The visibilities measured on the longest baseline (top) and closure phases (bottom) are plotted as a function of time over the night. For 
each block of five exposures (red data points), we compute the average observables and correct the calibrator visibilities to correct for their 
intrinsic UD diameters (black crosses). After rejecting data sets with strong intrinsic scatter, we average the calibrator measurements to compute 
the instrumental transfer function (grey line). 



signal-to-noise (SNR) ratio. Furthermore, we reject interfero- 
grams which were taken with an optical path delay (OPD) larger 
than 4 /mi in order to avoid systematic degenerative effects. To 
determine the transfer function, we first plot the observables for 
a representative spectral channel as a function of time (Fig. 2). 
Then, we reject observation blocks which show a significant in- 
trinsic scatter between the five individual exposures, typically 
indicating some degraded quality due to poor atmospheric con- 
ditions. While very few observations have to be rejected for the 
K-band, only a few if -band measurements provide reliable re- 
sults, reflecting the lower brightness of the object in this spec- 
tral window and the lower stability of the atmosphere at shorter 
wavelengths. Finally, we correct for the intrinsic diameters of the 
calibrator stars used (Tab. 1) and average the diameter-corrected 
calibrator observations for each night and for each spectral chan- 
nel in order to yield the transfer function. The calibration error 
is estimated from the scatter of the individual calibrator observa- 
tions over the night. As can be seen in Fig. 2, the calibration un- 
certainty is typically about ~ 3% for visibility and ~ 5° for CPs. 
Besides this calibration error, we include the statistical error es- 
timate provided by amdlib. Since this procedure results in un- 
realistically small error bars for some individual measurements, 
we add a constant calibration uncertainty of 3% for the visibility 
measurements. As can be seen in Fig. 2, the measurements ob- 



tained with different DIT values agree very well with each other, 
indicating good data quality. 

For some of our observations on RCrA, we employed the 
AMBER beam commutation device (BCD, Petrov et al. 2007). 
This calibration device allows one to exchange the beams of two 
of the three telescopes within a few seconds and to trace poten- 
tial drifts in the CP transfer function on much shorter time scales 
than a typical VLTI target/object cycle (~ 30 min). Furthermore, 
the BCD device is located close to the beginning of the optical 
path of the AMBER instrument and inverts the phase sign of 
the two exchanged beams, while the instrument-internal phases 
are not affected (Millour et al. 2008). Due to this optically intro- 
duced change in the CP sign, it is possible to distinguish real as- 
trophysical CP signatures from potential systematic instrumen- 
tal artifacts. Our BCD observations on RCrA were performed 
on 2008-06-02 and show the expected change of sign (Fig. 3), 
confirming the astrophysical origin of the detected CP signals. 

Fig. 4 shows the position angle dependence of the derived 
visibilities and CPs for two spectral windows around 2.0 and 
2.2 //rn. In order to interpret these position angle dependent vari- 
ations, it is necessary to model the object morphology, taking 
also the variations of the projected baseline length with position 
angle into account, as will be done in Sect. 4. In Fig. 5 (left), 
we plot the visibilities measured for the PA range 83-113° as 
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Fig. 5. Left: Visibilities derived from four RCrA observations taken on 2008-06-09 and covering the PA range between 83° and 113°. For the 
H- and /T-band spectral channels covered by our observations, the measured visibilities follow practically an identical visibility profile. Right: 
Comparing the visibility measured during four observing nights towards very different position angles (coded with different colors) reveals a 
strong position angle dependence of the visibility function, in particular at spatial frequencies above 8 X 10 6 i. 
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Fig. 6. Left: Comparison of the /T-band visibilities measured for the PA range 83-1 13° with various point- symmetric model geometries. It is evident 
that one-component models, such as UD, RING, and GAUSS geometries, are not consistent with the visibility profile, while a point- symmetric two- 
component Gaussian model (2-GAUSS) can reproduce the visibilities reasonably well, possibly indicating the presence of two spatial components 
such as an envelope and a disk (Sect. 4.2). Right: Fitting the 2-GAUSS model to the visibility data measured towards very different PAs shows that 
this model cannot reproduce the measured position angle dependence of the visibilities, in particular on the longer baselines. 



a function of spatial frequency. The spatial frequency v = B/A 
provides a measure for the resolving power achieved with a cer- 
tain observation and is therefore proportional to the projected 
baseline length and inverse proportional to the observing wave- 
length A. As can be seen in the left panel, the visibilities mea- 
sured for different spectral channels in the H- and /T-band can 
be very well represented with the same visibility profile. This 
indicates that the brightness distribution does not show a strong 
wavelength dependence 2 , which has interesting implications on 
the temperature-distribution of the emitting physical structure 
(see discussion in Sect. 5.3). Given the relatively small number 
of reliable //-band measurements, we consider only the ZT-band 
measurements for the following quantitative modeling. Also, the 
measured //-band closure phases are associated with very large 
error bars and are therefore omitted for the quantitative analysis. 



2 This argument is valid only for the visibility function V, but might 
not apply to the wavelength dependence of the closure phases, which are 
more sensitive to small-order changes in the object morphology (V oc 
v 2 , O oc y 3 ; see Lachaume 2003). 



4. Modeling 

In this section we describe the geometric and physical mod- 
els which we fitted to our interferometric data to constrain 
the spatial distribution of the circumstellar material around 
RCrA. Besides circumstellar dust and gas emission, the K- 
band might contain flux contributions from the stellar photo- 
sphere. Therefore we estimate the photospheric flux contribu- 
tions by comparing the de-reddened SED of R CrA with atmo- 
spheric models (Fig. 12), yielding relatively small values be- 
tween ~ 0.5% (assuming the spectral type F5, Garcia Lopez 
et al. 2006) and ~ 3% (spectral type B8, Bibo et al. 1992) of 
the total flux for the ^-band. Given this result, it seems justi- 
fied to neglect the photospheric contributions and to include only 
the thermal emission of hot circumstellar material in our mod- 
eling process (with the exception of the curved rim model, see 
Sect. 4.6). 

In order to fit the described models to our interferomet- 
ric data, we first generate model images for each model and 
each set of parameters. By construction, the geometric models 
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Table 2. Model-fitting results. 
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132 16 


0.60 


0.35 




27 


0.68 


2.30 


6.16 


3.25 


CURVED RIM 






180 35 








32 


0.50 


1.77 


3.28 


2.14 



Notes - (a) This column gives the model orientation, measured East of North. For the BINARY model, gives the PA of the separation 
vector, while for the SKEWED RING, VERTICAL RIM, and CURVED RIM models, the orientation of the ellipse/rim major axis is given. Due 
to a lack of closure phase calibration observations on the E0-G0-H0 array configuration, we are unfortunately not able to unambiguously define 
the CP sign, resulting in a 180° -ambiguity in the derived position angles, (b) In the fitting process, this parameter was fixed, (c) This model is 
point- symmetric, resulting in a CP which is identical zero. 
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Fig. 4. Visibilities (top) and closure phases (bottom) measured for 
RCrA with the E0-G0-H0 telescope configuration towards different 
position angles and for two representative spectral bands around A = 
2.0 11m (1.96 - 2.04 11m) and 2.2 jjm (2.16 - 2.24 11m). In both cases, 
three spectral channels around the central wavelength were averaged. 



(Sects. 4.1 to 4.5) are monochromatic, while for our curved rim 
model (Sect. 4.6), we compute the brightness distribution for 
each spectral channel separately, assuming the computed surface 
layer temperature for each disk annulus. From these model im- 
ages (with a scale of 0.2 mas/pixel), we compute the Fourier am- 
plitudes (visibilities) and Fourier phases for the wv-coordinates 
covered by the interferometric observations. Model visibilities 



Star 


V 


K 


Spectral 


dun 










Type 


[mas] 




HD 108570 


6.14 


4.12 


K1III 


0.72 + 0.05 


(a) 


HD 101328 


7.44 


3.75 


K4III 


1.00 + 0.01 


(b) 


HD 104479 


4.75 


3.92 


K0III 


0.76 + 0.01 


(a) 


HD 106248 


6.35 


3.73 


K2III 


0.94 + 0.01 


la) 


HD 111123 


1.30 


1.98 


B0.5IV 


0.84 + 0.06 


(a) 


HD 121384 


6.01 


4.00 


G8V 


0.76 + 0.05 


(a) 


HD 135452 


6.88 


3.43 


K3III 


1.13 + 0.02 


(a) 


HD 145921 


6.15 


3.62 


K2III 


0.96 + 0.01 


ib) 


HD 152534 


6.77 


4.54 


G8III 


0.60 + 0.04 


(a) 


HD 170499 


7.76 


3.05 


K4III 


1.24 + 0.02 


ib) 


HD 174631 


6.11 


2.94 


K1III 


1.25 + 0.02 


ib) 


HD 177756 


3.43 


3.56 


B9V 


0.53 + 0.04 


(a) 


HD 181110 


7.30 


3.71 


K3III 


0.95 + 0.01 


ib) 


HD 183925 


6.75 


2.92 


K5III 


1.44 + 0.02 


(a) 



Notes - The V-band magnitudes were taken from SIMBAD and 
the ^T-band magnitudes from the 2MASS point source catalog. 

References - (a) UD diameter computed with ASPRO. (b) UD 
diameter determined with the ESO CalVin tool. 



V and closure phases O r are computed for each spectral chan- 
nel separately and then compared to the observables (V, O) and 
their uncertainties (cry, cr®) by adopting^ = x^y + Afro as me 
likelihood estimator, with 



X 



^ V " Ny ^ 



= 1..JV 



1 - — y 



= 1..JV 




(1) 

(2) 



where N v and are the degrees of freedom for the individual 
visibility and closure phase measurements, respectively. In order 
to find the best-fit set of parameters for each model, we vary all 
free parameters on a parameter grid and search for the global x 2 r 
minimum. 
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In the following, we give a detailed description of the ap- 
plied models and show the obtained best-fit results for the spec- 
tral channel around 2.2 pm (Figs. 5 to 12). A comparison be- 
tween the model and the data obtained in all spectral channels 
is shown in Figs. 13 and 14. The derived best-fit parameters are 
given in Tab. 2. 

4.1. UD, RING, GAUSS: Point-symmetric models 

For a first quantitative interpretation of our data, we employed 
the commonly used uniform disk (UD), ring (RING), and 
Gaussian (GAUSS) geometries. By construction, these simple 
geometries are point-symmetric and, thus, not able to reproduce 
any object elongation or object asymmetries. Therefore, in con- 
trast to the following Sects. 4.3 to 4.6, here we aim only at re- 
producing the measured visibilities, being aware that the real 
object geometry is, in fact, not point- symmetric. Nevertheless, 
these simple models can be useful to estimate the characteris- 
tic size of the source brightness distribution and to allow a 
comparison with interferometric studies of other objects. In or- 
der to minimize potential influences from object elongation, we 
applied these model fits to a data subset, covering a smaller PA 
range. Due to the large number of available independent mea- 
surements, we have selected the PA interval 83-113°. From the 
obtained best-fits (grey curves in Fig. 6, left) it is evident that 
neither of these point- symmetric standard geometries can repro- 
duce the measured visibilities. 

4.2. 2-GAUSS: Indications for a disk+envelope geometry 

A plausible explanation for the detected strong deviations be- 
tween standard model geometries and the measured visibilities 
might be the presence of multiple components, such as a disk 
and an envelope component. This scenario also gains support 
from earlier studies which required the presence of a particularly 
massive circumstellar envelope (Natta et al. 1993) to reproduce 
the SED of RCrA. 

As a simple two-component geometry, we considered a 
model which consists of two Gaussian components. In a first 
attempt (Model 2-GAUSS), we fixed the center of the two 
Gaussians and varied only their FWHM diameter (0\, 62) and 
intensity ratio (h/h), which allows us to reproduce the visi- 
bility profile measured towards a narrow PA-range reasonably 
well (Fig. 6, left). The parameters of this two-Gaussian model 
(model 2-GAUSS) are listed in Tab. 2. Of course, being point- 
symmetric, the model cannot reproduce the measured non-zero 
closure phases. However, the model allows one to explain the 
pronounced change of slope, which can be seen in the visibility 
function around spatial frequency v ~ 8x 10 6 A. At shorter spatial 
frequencies, the measured visibilities do not show any significant 
position angle dependence, while at higher spatial frequencies, 
the visibilities differ significantly toward different position an- 
gles (see Fig. 5, right). This strongly suggests the presence of 
at least two spatial components, one being rather extended and 
point- symmetric (dominating at v < 8 x 10 6 A) and the other be- 
ing more compact and strongly asymmetric (and dominating at 
y > 8 x 10 6 /i). Based on this rather general argument, we include 
two spatial components in all model geometries, namely 

(a) an extended, symmetric component, which we represent 
with a Gaussian (in the following, this component will be 
referred to as "Envelope") and 

(b) a compact, elongated, and asymmetric component (in the fol- 
lowing referred to as "Disk"), which has the purpose of re- 
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Fig. 7. Simulation of the expected visibility (top) and closure phase 
(bottom) signatures of a hypothetical companion star, plotted as func- 
tion of the binary flux ratio. For these simulations, we assume a point- 
symmetric extended component (for which we use the 2-GAUSS ge- 
ometry discussed in Sect. 4.2) and add a compact companion at a sepa- 
ration of 60 mas, corresponding to the minimum separation of the pro- 
posed wide- separation companion star (Takami et al. 2003). Then, we 
compute the residuals between the model with and without companion 
star and measure the amplitude of the wavelength-differential visibil- 
ity and CP modulation for the ^T-band spectral window. The dashed 
horizontal lines give the achieved wavelength-differential visibility and 
closure phase accuracy, indicating that our observations should be sen- 
sitive to any companion star contributing more than q = 1 : 40 of the 
total flux (shaded area). 



producing the measured non-zero closure phases and the de- 
tected position angle dependence of the interferometric ob- 
servables. 

Of course, for our modeling, the two components only represent 
geometric quantities on different spatial scales, and their deno- 
tion as "Envelope" and "Disk" is, at the modeling stage, done 
only for the sake of clarity. Therefore, we include in all follow- 
ing models (with the exception of the 2-GAUSS and BINARY 
model, which consist of two Gaussian components itself) an en- 
velope component, which is represented by two free parameters, 
namely the FWHM diameter (0 2 ) and the envelope/disk flux ra- 
tio(/ 2 //i). 



4.3. BINARY: Asymmetric two-component Gaussian model 

The presence of close companions can introduce strong asym- 
metries in the source brightness distribution and might cause 
strong non-zero closure phases such as those detected in our 
VLTI/ AMBER measurements. For RCrA, this interpretation 





Fig. 8. Bispectrum speckle image reconstructed from //-band speckle 
interferograms recorded with the ES0 3.6m telescope on 2007-04-24 
(field-of-view 6"). At the achieved diffraction-limited angular resolu- 
tion of 95 mas, RCrA appears as a point-source. The first and second 
diffraction ring are clearly visible with noise-features on the 1 %-level, 
allowing us to rule out the existence of a companion star in this separa- 
tion/brightness contrast range. 



seems particularly appealing, since two earlier studies proposed 
the presence of multiple stellar sources to explain the mea- 
sured spectro-polarimetric signatures (Takami et al. 2003) and 
the unexpected hard X-ray emission (Forbrich et al. 2006). From 
their spectro-astrometric signal, Takami et al. (2003) derived a 
lower limit of 8 AU (> 60 mas) for the apparent separation of 
this hypothetical companion. Given this wide separation, the bi- 
nary should be clearly detectable in our AMBER data. As il- 
lustrated in Fig. 3 of Kraus et al. (2009), binaries with such a 
wide separation cause a high-frequency cosine modulation in the 
wavelength-differential visibilities and closure phases (this argu- 
ment holds as long as the separation is large compared to the di- 
ameter of the individual components). Neither the measured vis- 
ibilities nor closure phases show such a systematic, wavelength- 
differential modulation (Fig. 13 and 14), which clearly indicates 
that no companion with a separation p > 20 mas is significantly 
contributing to the near-infrared emission. In order to investigate 
up to which flux ratio q = / C ompMotai our observations would 
be sensitive to the presence of such a hypothetical compan- 
ion star, we simulated the K-band wavelength-differential sig- 
natures of a compact companion star around a point- symmetric, 
extended component such as derived in Sect. 4.2 and compare 
the amplitude of the predicted wavelength-differential modula- 
tions with the achieved differential visibility and CP accuracy. 
Based on these simulations (Fig. 7), we estimate that our ob- 
servations rule out the existence of a companion star for the 
separation range 20 mas < p < 200 mas and for ^T-band flux 
ratios brighter than -1:40. The upper limit in this separation 
range is given by the field-of-view of the used VLTI/AT tele- 
scopes. Besides these constraints from VLTI/AMBER interfer- 
ometry, we obtained bispectrum- speckle interferometric obser- 
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vations with the ESO 3.6m telescope (//-band; data from 2007- 
04-24 and 2007-04-27) using the Rockwell HAWAII detector of 
our visitor speckle camera. For image reconstruction we used 
the bispectrum speckle interferometry method (Weigelt 1977, 
Weigelt & Wirnitzer 1983, Lohmann et al. 1983). Both the de- 
rived power- spectrum and the reconstructed diffraction-limited 
bispectrum speckle image (Fig. 8), shows RCrA as a point- 
source, ruling out the existence of a companion star at separa- 
tions p > 60 mas, down to a flux ratio of 1:40. 

In order to test for the presence of a companion at shorter 
separation, we also fitted a binary star model to our AMBER 
data. Since it is likely that the near-infrared emission of the two 
components is not dominated by photospheric emission, but in- 
stead by the thermal emission of hot circumstellar material, we 
represent the two components not with point-sources, but with 
extended geometries, namely Gaussians. The choice of this ge- 
ometry is motivated by our results from Sect. 4.2, where we 
showed that two Gaussians are well suited to represent the ra- 
dial intensity profile towards some PA ranges. 

In this model, the free parameters are the FWHM sizes of 
the two Gaussians (6\, 62), their flux ratio (/2//1), their angular 
separation (p), and the position angle (0). As shown in Fig. 9, 
the best-fit binary model can provide only a moderate fit to the 
measured visibilities and CPs (xj = 3.93). Furthermore, as will 
be discussed in Sect. 5.2, the found best-fit solution is not phys- 
ically meaningful in the context of a companion star scenario. 

4.4. SKEWED RING model 

For passive irradiated disks, most of the near-infrared emission 
should emerge from hot dust located in a narrow annulus around 
the dust sublimation radius. Therefore, Monnier & Millan-Gabet 
(2002) and others have argued that ring geometries might pro- 
vide a good approximation for the appearance of YSO disks on 
AU-scales. While for face-on disks, point- symmetric rings seem 
appropriate, ring models are likely to be an over-simplification 
for disks seen under significant inclination. Since one side of the 
rim should appear brighter than the other, Monnier et al. (2006) 
proposed a simple mathematical modification of the standard 
ring model using a sinusoidal modulation of the ring brightness 
as a function of the azimuthal angle. We apply a similar model- 
ing approach using the intensity distribution 

/ 4ln2 - (r - pR) 2 \ 
/(r, a) = (1 - s • sin(a - 0)) • exp I K p 1 , (3) 

where the first term modulates the brightness distribution as a 
Gaussian centered at radius R and of fraction width / (FWHM). 
The brightness distribution /(r, a) is defined in a polar coordi- 
nate system, where r denotes the radius from the center and a 
the azimuth angle. The projection factor p can be used to take 
inclination effects into account and is given by 

p = cos / (cos 2 / • cos 2 (a - 0) + sin 2 (a - 0)) / . (4) 

In both equations, denotes the PA along which the major axis 
of the apparent disk ellipse is aligned. The fitted free parameters 
are the ring radius (R), the fractional ring width /, the position 
angle 0, the inclination /, the skew parameter s, as well as the two 
envelope parameters (62, hlh) described above. For this skewed 
ring model, our fitting algorithm can reach good agreement with 
the data (x 2 r = 1.96; Fig. 10). 



4.5. VERTICAL RIM model 

Over the last decade, detailed physical models have been devel- 
oped to describe the 3-D structure of the inner dust rim around 
Herbig Ae stars. Dullemond et al. (2001) provided a first math- 
ematical description of the inner rim geometry and considered 
a perfectly vertical inner rim. Since a vertical rim should result 
in very strong asymmetries, it seems promising to test whether 
such a model can explain the strong closure phases detected in 
our R CrA data. To compute the brightness distribution, we as- 
sume that all near-infrared flux emerges from the hot, illumi- 
nated rim surface. In this model, the free parameters are the rim 
radius R, the rim scale height H/R, the position angle 0, the sys- 
tem inclination /, and the envelope parameters 62 and hlh- Our 
resulting best fit is shown in Fig. 11, but provides no satisfactory 
representation of the data (x 2 r = 3.25). 

4.6. CURVED RIM model 

Following the pioneering work by Natta et al. (2001) and 
Dullemond et al. (2001), several studies aimed to refine the 
theoretical description of the rim geometry. For instance, IN05 
pointed out that the dust sublimation temperature should depend 
on the local gas density. Accordingly, the dust sublimation tem- 
perature should be highest in the disk midplane and decrease 
with scale height //, resulting in a curved shape of the inner rim. 
In order to further increase the physical accuracy, Tannirkulam 
et al. (2007) included the effect of dust grain sedimentation, 
showing that the presence of a population of small and large dust 
grains results in an even stronger rim curvature than predicted by 
the IN05 model. 

To compute the radial dependence of the scale height of the 
disk photospheric layer H/R, we follow the analytical approach 
by IN05. The dependence of the dust sublimation temperature 
on the local gas density p g (in g/cm 3 ) is described by the relation 
from Pollack et al. (1994), namely r sub i = 2000 K • p° 0195 . In 
order to avoid unrealistic cutoffs at the outer edge of the puffed- 
up inner rim, we include the more extended disk regions in our 
model as well and construct our disk with the following compo- 
nents: 

(a) Puffed-up inner rim: This region extends from R^\ to the 
point where the rim surface layer becomes optically thin to 
the stellar radiation. 

(b) Shadowed region: Entering the optically thin regime, the rim 
can no longer maintain its puffed-up shape and the scale 
height decreases as d(H/R)/dR = -1/(87?), as derived by 
Dullemond et al. (2001). 

(c) Flared outer disk: At larger radii, the disk might en- 
ter a flared shape, which we parameterize with H/R = 
H (R/R ) 9/S and T(R) oc 7T 3/4 (Kenyon & Hartmann 1987). 

In addition, the model includes the photospheric emission, which 
we estimate from the SED by comparing the stellar flux with the 
measured near-infrared flux. Please note that our near-infrared 
observations are mainly sensitive to regions (a) and (b), while the 
parameters for the outer flared disk region are not constrained by 
our model fits. Therefore, we fixed the parameters H Q and Rq to 
the typical values of // = 0.1 AU and Ro = 10 AU. As total disk 
mass, we assume M d i sk = 0.012 Mq (Groppi et al. 2007) with 
a standard gas-to-dust ratio of 100:1. Besides the luminosity L 
which is incident on the disk rim, the most sensitive parame- 
ter in the IN05 model is e. This variable is defined as the ratio 
of the Planck mean opacity at the dust evaporation temperature 
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Fig. 10. Left: Comparison between the AMBER observables and the model visibilities (top) and closure phases (bottom) from our SKEWED 
RING model. Right: Image corresponding to our best-fit model. 



^subi to the Planck mean opacity at the stellar effective temper- 
ature T+ and defines the cooling efficiency of the dust grains. It 
therefore depends on the precise dust chemistry and grain size 
distribution. For a given grain population, e typically increases 
for larger grain sizes, moving the rim location closer to the star 
(see IN05). Above a certain threshold (e > e cr = 1/ V3), the lo- 
cation of the optically thick dust rim stays nearly constant, and 
the largest grains can only survive closer to the star in a opti- 
cally thin regime. Therefore, we imposed for our fit that e cannot 
exceed e cr . In addition, the model assumes that absorption and 
continuum emission by gas inside the dust sublimation radius is 
negligible. This assumption is supported by the theoretical work 
of Muzerolle et al. (2004), who found that for the accretion rate 
of RCrA (M acc = 10" 712 M Q /yr, Garcia Lopez et al. 2006), the 
gas located inside the dust sublimation radius should be optically 
thin. 

After computing the disk scale height H (r) and blackbody 
temperature T(r) for each disk annulus, we construct synthetic 
images of the disk brightness distribution. By computing these 
images for each spectral channel separately, we can simulate the 
wavelength-differential changes resulting from the disk temper- 
ature distribution. 

The free disk model fitting parameters are the position angle 
0, the system inclination i, the grain cooling efficiency 6, and 
the envelope parameters 62 and hlh - Of course, another crucial 
parameter is the luminosity L, which is heating the inner dust 
rim, and which might contain stellar light contributions (L*) as 
well as contributions due to active accretion (L acc ). Although the 
precise amount of either of these contributions is still not well 
known, it is likely that the total incident luminosity L := L*+L acc 
is neither as low as ~ 2 Lq (corresponding to the stellar lumi- 



nosity of a F5 star), nor that the full bolometric luminosity of 
~ 99 L© (Bibo et al. 1992) derived from KAO photometry can 
be solely attributed to this single source. Therefore, we treat L 
as a free parameter and keep the other stellar parameters, which 
have a less significant influence on our modeling results, fixed 
to the values corresponding to a B8-type star (T^ = 12, 000 K, 
Af* = 3.0 M Q , Bibo et al. 1992). Then, we adjust the absorp- 
tion coefficient A v to roughly reproduce the SED of R CrA at 
UV/visual wavelengths. Following Bibo et al. (1992), the total- 
to- selective extinction ratio Ry was fixed to 4.7. In order to in- 
clude the emission of the puffed-up inner rim in our model SED, 
we follow the approach outlined by IN05, assuming that the rim 
surface layer reprocesses the incident light and re-emits it as a 
blackbody of temperature ~ 1400 K. When modeling the SED 
of RCrA, one faces the general problem that the visual/near- 
infrared emission also contains major scattered flux contribu- 
tions from the circumstellar envelope, as described in Sects. 4.2 
and 5.1. Due to this strong "contamination", we decided not to 
include the SED as an additional modeling constraint, but to 
treat the fraction of stellar luminosity which is reprocessed at 
near-infrared wavelengths := Lnir/L as a free model parame- 
ter, which is adjusted to match the measured SED of R CrA. In 
Sect. 5.3, we will compare the derived value with the predic- 
tions from the IN05 model and discuss the consistency between 
the measured SED and the puffed-up inner rim scenario. 



Using our modeling approach, we find that the parameter 
combination L = 29 Lq and e = e CY is able to simultaneously re- 
produce the measured visibilities and closure phases (x 2 r = 2.14, 
Fig. 12, top, left) as well as the visual- to near-infrared SED 
(with Ay = 5.0 and f3 - 0.35; see Fig. 12, bottom, left). For 
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Fig. 11. Left: Comparison between the AMBER observables and the model visibilities (top) and closure phases (bottom) from our VERTICAL 
RIM model. Right: Image corresponding to our best-fit model. 



the following reasons, both parameters appear to be relatively 
well constrained: 

- Assuming a higher value for L would push the dust sublima- 
tion radius outwards, resulting in too low visibilities (e.g., for 
L = 50 L©, the best-fit^ -value increases already to 2.85). 

- A lower luminosity would require to increase J3 above phys- 
ically reasonable values (e.g., for L = 20 L©, already more 
than 60% of the stellar light would have to be reprocessed by 
the inner rim). 

- 6-values significantly lower than 6 cr seem to be excluded, 
since this would again result either in too low visibilities, 
or require to increase J3 substantially. 

We conclude, that already with a relatively small number of free 
parameters (L, e, 0, i, 62, h/h), it is possible to simultaneously 
reproduce the interferometric observables (x 2 r = 2.14) and the 
relevant parts of the SED. At wavelengths > 3 yum, our model 
is not able to reproduce the SED, since we do not include the 
flux contributions from the more extended disk regions nor from 
the circumstellar envelope. In Fig. 12 we show the model vis- 
ibilities and closure phases, the SED, as well as the computed 
photospheric scale height H/R and the temperature distribution 
for our best-fit model. In this model, the inner rim radius is lo- 
cated at 0.38 AU. 



5. Discussion 

5.1. Radial intensity profile 

As shown in Sect. 4. 1 , it is evident that the radial intensity profile 
measured toward R CrA cannot be well represented with simple 



one-component geometries, such as rings. Instead, the visibility 
function has an approximately linear slope (Fig. 5, left), showing 
remarkable similarities with the visibility profile measured by 
Acke et al. (2008) on the much more massive YSO MWC 297. 
In both cases, the visibilities could be approximated reasonably 
well with two Gaussian components. For the case of R CrA, we 
interpret this finding with the presence of an extended optically 
thin envelope, which contributes about one-third of the total K- 
band flux, likely through scattered light. This extended emission 
might represent scattered light from an optically thin spherical 
halo or from the walls of an associated outflow cavity. The pres- 
ence of an envelope around R CrA was already deduced earlier 
from SED fits (Natta et al. 1993) and from polarimetric observa- 
tions (Clark et al. 2000) and seems in line with the early evolu- 
tionary stage of the object. 



Independent of the precise physical interpretation of the 
measured visibility profile, our finding shows that the commonly 
applied ring model fits often represent an over-simplication of 
the complex environment around YSOs. Since these ring fits 
are commonly applied to single-baseline, broad-band interfer- 
ometric observations, the uncertainties of the real underlying ra- 
dial intensity profile likely contribute significantly to the scatter 
which is observed in size-luminosity diagrams (Monnier et al. 
2005). Therefore, very detailed interferometric studies on some 
additional sources will be essential to identify the real underly- 
ing source structure and will have a direct impact on the proper 
interpretation of the data obtained in survey-type observations. 
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Fig. 12. Top, left: Comparison between the AMBER observables and the 
model visibilities and closure phases from our CURVED RIM model. 
Top, right: Image corresponding to our best-fit model for the wave- 
length 2.2 yi/m. Middle: Photospheric height (left) and temperature dis- 
tribution (right) of the inner rim for our best-fit model. Bottom, left: 
SED of R CrA including photometric data from Hillenbrand et al. (1992, 
data points) and an ISO spectrum (Frieswijk et al. 2007, SWS AOT-1 
HPDP, dark-red curve). The grey curve shows the photospheric emis- 
sion corresponding to a B8 model atmosphere (Kurucz 1970) with lumi- 
nosity 29 L©, while the dashed blue curve shows the reprocessed light 
from the rim. The black curve represents the sum of both components. 
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Fig. 13. Comparison of the measured wavelength-dependent AMBER visibilities and the model visibilities corresponding to the best-fit model 
assuming our BINARY model (Sect. 4.3), SKEWED RING model (Sect. 4.4), VERTICAL RIM model (Sect. 4.5), and CURVED RIM model 
(Sect. 4.6). 



5.2. Non-detection of the proposed binary companion 

Given the earlier speculations about the existence of a compan- 
ion star for RCrA (Takami et al. 2003; Forbrich et al. 2006, see 
discussion in Sect. 2), it is an important result that we did not 
find indications for binarity of RCrA from our study. As dis- 
cussed in Sect. 4.3, it is possible to rule out the presence of wide 
separation binaries (p > 40 mas) with rather general arguments. 
To probe for binaries with smaller separation, we constructed 
geometric models, which also take into account that each star 
might be associated with a circumstellar disk (represented by 



Gaussians in our model). Since this model provided only a poor 
representation of our data, we consider such a close binary sce- 
nario rather unlikely. Furthermore, the best-fit parameters of the 
binary model do not correspond to a physical solution, since the 
extension of the Gaussian components (1.9 and 0.8 AU) exceeds 
even their projected separation (0.8 AU), resulting in an overlap 
of the two components. Assuming that the physical separation 
is of similar order as the projected separation, a hypothetical bi- 
nary on such small spatial scales would not be stable and would 
quickly disrupt the circumstellar disks around both components. 
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Fig. 14. Comparison of the measured wavelength-dependent closure phases 
SKEWED RING model (Sect. 4.4), VERTICAL RIM model (Sect. 4.5), and 
of the three baselines and the projected length of the longest baseline #eo-ho 
array, the other two baselines have projected length of #eo-go = #eo-ho/3 and 



BINARY model 
SKEWED RING model 
VERTICAL RIM model 
CURVED RIM model 



with the predictions from our best-fit BINARY model (Sect. 4.3), 
CURVED RIM model (Sect. 4.6). In each panel, the position angle 
is given. Given the equal spacing of the employed linear telescope 

#G0-H0 = 25 E o-Ho/3. 



Of course, with the existing data, it is difficult to rule out sce- 
narios which might involve three or more spatial components, 
such as, for instance, a binary system with circumstellar and cir- 
cumbinary disks. Also, our observations do not rule out the ex- 
istence of a very deeply embedded companion, which might not 
contribute significantly to the near-infrared emission. However, 
this scenario is rather unlikely, since the SED of R CrA peaks in 
the ^T-band (Fig. 12, bottom, left), which makes the existence of 
an additional strong, far-infrared-emitting component unlikely. 



Our non-detection of a binary component around R CrA sug- 
gests that the detected hard X-ray emission is not associated with 
a hypothetical Class I source (Forbrich et al. 2006), but with a 
single Herbig Ae star. This finding on R CrA is in line with the 
study by Hamidouche et al. (2008), who concluded from statis- 
tical arguments that the X-ray emissions of Herbig Ae/Be stars 
have an intrisic origin. Since neither coronal nor shock-excited 
X-ray emission is expected for A-type stars, more theoretical 
work is clearly required in order to identify the X-ray-emitting 
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mechanism of single intermediate-mass pre-main-sequence stars 
such as R CrA. 

5.3. Constraints on the rim geometry 

Given that our best-fitting results were obtained with the skewed 
ring model (Sect. 4.4) and the curved rim model (Sect. 4.6), we 
consider it most likely that our observations are directly trac- 
ing the asymmetries introduced by material at the inner trunca- 
tion radius of the dust disk around RCrA. This interpretation 
is also supported by the fact that the H- and ^-band visibilities 
follow within the error bars the same visibility profile (Sect. 3 
and Fig. 5, left), providing important information about the tem- 
perature distribution of the near-infrared emitting disk mate- 
rial. Earlier spectro-interferometric observations on two other 
Herbig Ae/Be stars (Kraus et al. 2008; Isella et al. 2008) mea- 
sured in the //-band a significantly higher visibility than in the 
ZT-band, which was interpreted as the presence of an optically 
thick, hot gaseous disk component located inside of the dust sub- 
limation radius. For RCrA, we do not measure this visibility in- 
crease towards shorter wavelengths, indicating that the disk ma- 
terial in the probed temperature range (> 1000 K) is located at 
similar distances from the star, likely in a narrow region around 
the dust sublimation radius. 

Concerning the detected asymmetries, it is interesting to 
compare our results with earlier studies on HD 45677 (O < 27°, 
Monnier et al. 2006) and AB Aur (O < 4°, Millan-Gabet et al. 

2006) , which also detected asymmetries in the inner region 
around these Ae- and Be-type stars. Comparing HD 45677 and 
RCrA, it is remarkable that the closure phases from the two 
objects could be represented well with a skewed ring model. 
However, a major difference between both objects concerns the 
evolutionary stage: HD 45677 is likely an evolved B[e] star, 
while RCrA is a young, actively accreting Herbig Ae star. 
Compared to the Herbig Ae star AB Aur, a major difference con- 
cerns the physical nature of the asymmetry. For ABAur, the 
modeling of Millan-Gabet et al. showed that the asymmetries 
around this object originate from a localized region within the 
disk (maybe a hot accretion spot), while our data seems more 
consistent with a diffuse structure around the dust sublimation 
radius. Another Herbig Ae/Be star where strong asymmetries 
were detected is LkHodOl. Near-infrared aperture-masking ob- 
servations by Tuthill et al. (2001) revealed a skewed circular 
structure at a distance of ~ 3.4 AU around this rather massive 
(Af* « 5 - 10 Mq) and luminous (L* « 480 Lq) Herbig Be 
star. Remarkably, the structure in their observations shows some 
similarities with our best-fit skewed ring model, which might 
suggest that, contrary to earlier suggestions (Vinkovic & Jurkic 

2007) , Herbig Ae and Herbig Be have a common rim structure. 
To determine the precise rim shape, we have fitted three dif- 
ferent models, namely the skewed ring model (Sect. 4.4), the ver- 
tical rim model (Sect. 4.5), and the curved rim model (Sect. 4.6). 
Formally, the best agreement was found with the geometric 
skewed ring model (xj = 1.96), although it should be noted that 
the skewed ring model is a purely geometric model and has the 
largest number of adjustable free parameters (5 disk & 2 en- 
velope parameters). On the other hand, the curved rim model 
(Xr = 2.14) has less free parameters (3 disk, 2 envelope, and 
1 stellar parameter) and is based on detailed theories about disk 
structure and dust properties, allowing a direct physical interpre- 
tation of the fitted parameters. In this context, the dust cooling 
efficiency parameter e is of particular interest, since this param- 
eter is directly related to the dust properties in the inner disk 
regions. For instance, Isella et al. (2006) used the curved rim 



model to derive the maximum dust grain size for a sample of 
five Herbig Ae/Be stars, finding in all cases relatively large dust 
grains (> 0.2 jim). Assuming the same Silicate dust chemistry, 
we derive for RCrA dust grain sizes larger than 1.2 fim (corre- 
sponding to 6 > 6 cr ), indicating that the dust in the inner disk 
regions was already significantly processed by grain growth. In 
order to reproduce the interferometric data and the SED simulta- 
neously, we require that the incident luminosity, which is heating 
the inner dust rim, is around 29 Lq and that 35% of the stel- 
lar light is reprocessed at near-infrared wavelengths. This frac- 
tion is somewhat higher than the theoretical value of ~ 20% 
predicted by the IN05 model (see Fig. 4 in IN05). However, 
this might reflect the fact that a significant fraction of the near- 
infrared emission is likely not reprocessed light from the in- 
ner rim, but scattered light contributions from a circumstellar 
envelope (Sects. 4.2 and 5.1). Although it is currently unfortu- 
nately not possible to properly decompose the disk and envelope 
SED, it is likely that the ft- value has to be corrected by a similar 
factor as the determined envelope/disk ZT-band flux ratio of 1/3 
(Sect. 4.2), which would result in a good agreement with the ex- 
pected theoretical value. Alternative explanations for the derived 
high value of J3 would include an underestimation of the gas disk 
mass (Mdisk = 0.012 Mq) and the resulting underestimation of 
the Silicate sublimation temperature or the presence of highly re- 
fractory dust species which sublimate at higher temperature than 
Silicate grains. 

For the derived luminosity of ~ 29 Lq, we consider that 
the dominant contribution is photospheric emission, while the 
contributions from active accretion are likely less significant, 
as indicated by the observed low Bry-line luminosity (Garcia 
Lopez et al. 2006). Accordingly, only about one-third of the to- 
tal bolometric luminosity (~ 99 Lq, Bibo et al. 1992) can be 
attributed to R CrA, while the remaining fraction is likely due to 
contamination from other, more deeply embedded members of 
the Coronet cluster in the large (~45") field-of-view of the KAO 
satellite. 

For the disk inclination angle, we consider the value de- 
termined with the curved rim model (i = 35°) our most reli- 
able estimate 3 and note that the value is also in good agre- 
meent with the inclination angle of ~ 40°, which was derived 
by Clark et al. (2000) for the reflection nebula NGC 6729 using 
polarimetric observations. NGC 6729 also shows some remark- 
able fine- structure, in particular two bow shock-like features (la- 
beled E and F in Clark et al. 2000), which appear east of RCrA 
(Fig. 15, left and middle). As already suggested by Clark et al., 
these structures might be created by vigorous, periodic outflow 
activity from RCrA. The sub-AU scale disk structure resolved 
by our AMBER observations is oriented approximately perpen- 
dicular to these bow shocks (0 = 180, 190° for the CURVED 
RIM model and the SKEWED RING model, respectively), sug- 
gesting that R CrA is the driving engine which has created these 
bow shocks (Fig. 15, right). The derived disk orientation is also 
consistent with the orientation of the polarization disk reported 
by Ward-Thompson et al. (1985, = 189 ± 5°; see Fig. 15, mid- 
dle). The bright feature extending from RCrA to the south-east 
(labeled C in Clark et al. 2000) was identified by Clark et al. as 
scattered light from the walls of a parabolic outflow cavity. 



3 In particular, this intermediate inclination angle seems more reli- 
able than the very low inclination of 14° determined with the skewed 
ring model, reflecting the fact that in the skewed ring model the domi- 
nant free parameter is the skew parameter s, while the inclination i has 
only a minor effect on the model appearance, and therefore, this param- 
eter is only poorly constrained. 
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Fig. 15. At visual wavelengths the field around RCrA is dominated by the reflection nebula NGC 6729 (left panel; color composite: blue: 5-band, 
green: V-band, red: 7?-band; North is up and east is to the left). Within the reflection nebula, two bow shock- like structures appear (middle panel, 
dashed lines), suggesting an east-western outflow axis, which is roughly perpendicular to the polarization disk reported by Ward-Thompson et al. 
(1985) and the sub-AU disk resolved by our VLTI/AMBER observations (right panel; CURVED RIM model image; please note that there is still 
an 180° -ambiguity in the derived orientation, as discussed in Tab. 2). The B/V/R-band images were taken with the Wide Field Imager (WFI) at the 
2.2-m MPG/ESO telescope on La Silla (ESO Press Photo 25a-b/00; image courtesy: European Southern Observatory and F. Comeron). 



When judging the quality of the obtained fit, it is interesting 
to note that the IN05 fitting results improved slightly by tak- 
ing the expected disk temperature gradient effects into account 
(Xr = 2.14). For comparison, when we simulate a monochro- 
matic brightness distribution (A = 2.2 /mi), we yield x 2 r - 2.38. 
Nevertheless, with a^r-value of 2.14, the quality of the best fit 
IN05 model is still relatively poor, indicating that the rim geom- 
etry might still not be adequately represented by the employed 
rim model. Therefore, it will be important to investigate in future 
studies how the models can be modified to yield a better rep- 
resentation of the measured visibilities and closure phases. For 
instance, including dust sedimentation and grain growth effects 
should result in an even more curved rim shape and a broader 
near-infrared emitting rim region than predicted by the current 
model (Tannirkulam et al. 2007). Based on the good fitting re- 
sults obtained with a skewed ring model with a large fractional 
ring width of 0.8 (Sect. 4.4), we expect that such a modifica- 
tion might improve the visibility/phase fit. Possibly, an improve- 
ment could also be obtained by including the optically thin dust 
condensation zone which is expected to extend inwards of the 
dust rim, yielding a more diffuse rim structure (Kama et al. 
2009). Other very promising modifications concern the presence 
of multiple grain species, in particular highly refractory metal 
oxides such as iron or corundum, and a consistent treatment of 
backwarming effects. Kama et al. have shown that these effects 
can move the dust sublimation radius nearly by a factor of 2 
closer to the star. With our fitting procedure, this change would 
result in a higher derived stellar luminosity (i.e. more consistent 
with the measured bolometric luminosity) and a lower fraction 
of reprocessed light (i.e. closer to the theoretical expectation). 
Thus, we expect that these model refinements aim in the right 
direction in order to further improve the physical consistency of 
the puffed-up inner rim scenario for R CrA, although a detailed 
data modeling study will be required in order to test whether the 
model refinements also affect the rim morphology, as required to 
yield a better;^ 2 -fit. Finally, one should take into account that our 
observations probe relatively small physical scales in a highly 
dynamical and complex environment. Therefore, it cannot be 
excluded that our interferometric observations might also be af- 
fected by local inhomogenities or astrophysical processes, which 
are not yet completely included in the employed disk rim mod- 



els, such as active accretion or the influence of outflow launching 
on the disk structure. 



6. Conclusions 

We summarize our findings as follows: 

- Using 24 VLTI long-baseline interferometric measurements, 
we studied the near-infrared emission from the Herbig Ae 
star R CrA and could clearly resolve the inner circumstellar 
environment in the H- and AT-band. Besides this spectral cov- 
erage, our observations also provide a good sampling of the 
visibility function (0.8 < V < 0.1) towards a wide range of 
position angles (PA=16-113°). 

- Even when considering only a single position angle, the mea- 
sured visibility profile cannot be represented with the com- 
monly applied one-component models, but indicates a more 
complex object geometry, likely including a disk and an en- 
velope component. 

- The measured closure phase signals (up O ~ 40°) clearly 
indicate a strongly asymmetric brightness distribution. 

- We find that the detected asymmetries are likely not related 
to the presence of a companion star, but directly trace the ver- 
tical structure of the disk around the dust sublimation radius 
(^subi ~0.4 AU). When seen under intermediate inclination, 
the increased disk scale height causes obscuration and shad- 
owing effects, which result in the observed asymmetries. To 
constrain the precise rim geometry, we tested three geomet- 
ric and physical models, including the generic skewed ring 
model and rim models with a vertical and a curved rim shape. 
Clearly, models with a smooth brightness distribution (i.e. 
the skewed ring model and the curved rim model) provide 
a much better representation of our data than the sharp rim 
structure predicted by the vertical rim model. 

- Confronting our data with the detailed mathematical descrip- 
tion of the rim structure presented by IN05, we find that we 
can reasonably well reproduce the measured SED, as well 
as the visibilities and closure phases with a disk inclination 
angle of 35° and an incident luminosity of 29 L©. Given the 
absence of indications for strong active accretion, we con- 
sider that this luminosity is mainly of photospheric origin. 
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For the dust grain size, we find that the presence of relatively 
large Silicate dust grains (> 1.2 //m) is required to obtain 
agreement between the model and our data. The derived disk 
position angle of ~ 180 - 190° agrees well with the orien- 
tation of the polarization disk (6 = 189 ± 5°). Perpendicular 
to the disk axis, two bow shocks appear in the associated re- 
flection nebula NGC 6729, suggesting that the detected disk 
is driving an outflow, which has shaped the bow shock-like 
structures. 

Presenting one of the most comprehensive high-angular res- 
olution studies on the inner structure of YSO disks, our findings 
provide strong evidence for the existence of a curved, puffed- 
up inner dust rim in Herbig Ae stars. However, additional work 
is needed to further constrain the detailed rim geometry and 
to distinguish the rim morphology from local brightness inho- 
mogenities, which might, for instance, be caused by hot accre- 
tion spots, spiral disk density patterns, dumpiness, or other tran- 
sient phenomena. Given the complexity one likely faces on the 
probed sub-AU scales, this task will urgently require the model- 
independent aperture synthesis imaging capabilities which are 
just becoming available for various infrared interferometric fa- 
cilities. 
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